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1.  Introduction.  W  auuld  lihe  to  Introduce  tho  problem  with  en  example. 
How  should  we  design  an  experiment  to  compere  4  test  treatments  with  a 
control,  using  18  experimental  units?  As  a  statistical  question  we  will 
not  be  able  to  answer  it  unless  it  is  asked  in  a  more  precise  manner.  To 
begin  with  we  need  to  postulate  a  model  for  the  response  observed  upon 
application  of  a  treatment,  test  treatment  or  control,  to  an  experimental 
unit.  ihls  paperfwe-  shall  consider  three  possible  models: 


- -  - - - - -  - - - - 

V 

l)  O-way  elimination  of  heterogeneity  model  in  which  all  experimental  units 
are  homogeneous  before  application  of  treatments:  \ 


y±i  "  N  ♦  fci  ♦  *ij 


(1.1) 


T-)  1-way  elimination  of  heterogeneity  model  in  which  experimental  units  can  be 
divided  into  several  homogeneous  blocks:  ,4>. 


Tij  "  f  ♦  *i  ♦  ♦  «ij 


(1.2) 


2-way  elimination  of  heterogeneity  model  in  which  the  experimental  units 
can  be  conceptually  arranged  according  to  rows  and  columns: 


r iji  -  p  ♦  t±  ♦  8j  ♦  P}  ♦  tij^ 


(1.3) 


Now  we  can  be  more  precise  about  what  we  mean  by  comparing  test 
treatments  with  a  control.  Our  goal  is  to  estimate  the  magnitude  of  each 
(t£  -  t, ).  Assuming  that  the  error  component  e  in  the  model  is  homo- 
scedastlc,  the  method  of  least  squares  will  be  used  to  estimate  the 


vswstm 


3- 


contrasts  (t^  -  t, ) ;  this  happen*  to  ba  th*  boat  linear  unbiased  estimator 
(ti  -  t,).  In  aaalfning  the  treataenta  to  experimental  unite  we  have  to  make 
aura  that  the  contrasts  (t^  -  tt )  are  estimable.  In  eaae  we  have  more  than 
one  choice  for  making  this  assignment  ire  want  to  select  one  which  guaran¬ 
tees  high  efficiency  in  the  sense  of  achieving  th*  minimum  value  of 

4 

J  var  <ti  -  t.)  (1.4) 

or 

l!il4  <*1  "  °’5) 

A  design  which  gives  th*  minimum  in  (1.4)  will  be  called  an  A-optlmal 

design  and  one  which  gives  the  minimum  in  (1.5)  will  be  called  an  MV- 

optimal  design. 

Without  further  ado  we  give  deaigns  which  are  A-  and  MV-optlmal  under 
each  of  the  three  models: 

A-  amd  MV-optlmal  dealgm  under  model  (1.1): 

Assign  3  experimental  units  to  each  of  th*  4  test  treatments  and  6  to 
th*  control. 

A-  and  MV-optlmal  design  under  model  (1.2),  when  there  are  6  blocks  of  sis* 

3  each: 

Take  each  column  of  the  following  array  as  a  block: 

0  0  0  0  0  0 

1112  2  3 

2  3  4  3  4  4 


Her*  0  denotes  th*  control  and  1,  2,  3  and  4  the  test  treatments 


1  (1.3),  where  there  art  3  rows  and 


6  coluana: 

Aiiltn  tha  treataents  according  to  tha  following  array: 

1  0  3  4  2  0 
0  3  4  2  0  1 
4  2  0  0  1  3. 

This  aaall  axaaple  daaonatrataa  tha  fact  that  tha  problaa  of  finding 
an  optiaal  daaign  la  a  difficult  ono.  During  tha  paat  aavaral  yaara  thara 
haa  baan  a  concantratad  affort  to  idantify  and  conatruct  optiaal  daaigna 
for  tha  ganaral  problaa  of  coaparlng  v  teat  traataanta  with  a  control. 
While  in  thaaa  atatlatical  aattlnga  A*  and  MV-optlaality  are  tha  two  aoat 
natural  optiaality  criteria,  wa  aay  want  to  conaldar  other  optiaality 
criteria  aa  wall  in  certain  altuationa.  However,  in  our  view  tha  pub¬ 
lished  literature  on  othar  optiaality  criteria  baa  not  reached  a  level  of 
generality  for  aunaarlaatlon.  In  thla  paper  wa  shall  atteapt  to  suaaarise 
results  on  A-  and  MV-optlaal  designs,  which  wa  hope  will  be  useful  to  both 
tha  theoretician  and  tha  practitioner. 

In  Sections  2  and  3  wa  give  ganaral  results  for  A-  and  MV-optiaal 
designs  for  coaparlng  v  teat  traataanta  with  a  control  in  each  of  the  3 
aodala  (1.1),  (1.2)  and  (1.3).  In  Section  4  wa  give  nodal  robust  A-  and 
MV-optlaal  designs.  In  Section  5  we  suggest  various  approaches  for  finding 
efficient  designs  in  those  cases  where  A-  and  MV-optiaal  designs  are 
unknown.  In  Section  6  we  give  A-  and  MV-optiaal  designs  for  coaparlng  test 
treataents  with  two  or  aora  controls.  In  Section  7  we  give  an  overview  of 
the  literature  of  optiaal  designs  for  coaparlng  test  treataents  with 


controls. 
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2.  A«optlMl  Pt»ltn».  V*  shall  |l*i  A-optimal  daaigna  for  comparing  v 
teat  treatments  with  a  control  aaparataly  for  tha  O-way,  1-way,  and  2-way 
allalnation  of  hataroganalty.  Throughout  this  auction  tha  control  will  ba 
danoted  by  tha  syabol  0  and  tha  taat  traataanta  by  1,  2,  .  .  . ,  v. 


2.0.  A-optiaal  daaigna  for  0-way  alialnation  of  hataroganalty.  Our 
statistical  aat-up  consists  of  n  experimental  units,  and  our  aodal  of 
response  under  a  design  d  is 

ydlj  ■  H  ♦  *ij  (2.1) 

where  J  •  1,  .  .  .,  r^^,  i  ■  0,  1,  ....  v.  Hera  r,n  is  tha  number  of 
experimental  units  receiving  treatment  1.  He  assume  tha  model  to  be 
additive  and  homoscedastlc .  Tha  symbols  in  equation  (2.1)  have  their 
standard  meaning.  The  A-optlmal  design  minimizes 

Jj  <l/*ds  +  l'*di> 

subject  to  tha  restriction  rj,  +  rj,  ♦  .  .  r«jv  -  n.  Using 

elementary  mathematics,  it  is  easily  seen  that  an  A-optiaal  design  d 
will  have  Ir^*!  -  tj*j|  S  1  for  i,J  •  1,  .  .  .,v  and 

n  -  v(f,  ♦  $)  S  rj*0  in-  fa  —A  (2.2) 

^  J 

where 

*  H(2nv  ♦  v  -  v*)  -  (v*  ♦  v*  -  2v3  ♦  4n*v)i}/2v(v  -  1)M*  I  (2.3) 
Here  (x]  denotes  the  integer  part  of  the  decimal  expansion  for  x  >  0.  In 
the  case  v  is  a  square  and  n  »  m(v  +  vi)  for  an  integer  a,  the  A-optimal 
design  d*  is: 

rd*i  *  •  •  •  ■  rd*v  *  ».  rd*P  “  ■  vK 


Hi 


6 


2.1.  A-optlaa 1  d»>ltni  for  1-vay  elimination  of  h>t>ro«»n»lty.  Out 
statistical  sst-up  consists  of  b  blocks  of  size  k  each,  and  the  aodel  of 
response  under  a  design  d  1st 

Pdijp  -  P  ♦  ti  ♦  Bj  ♦  «ijp. 

where  i  •  0,  1,  .  .  v;  j  ■  1,  .  .  . ,  b  and  p  ■  0,  1,  .  .  n^j. 
Hera  n^^j  la  the  nuaber  of  tlaes  treatment  1  la  used  in  block  j.  Let  H4 
denote  the  aatrix  (n^j ) , 

b 

rdi  "  °dlj* 

C(|  ■  Dlag  ( rj, ,  rji ,  .  .  • ,  r^y)  “  k  *  , 


P  being  a  v  x  (v  ♦  1)  aatrix.  Then  an  A-optiaal  design  ainialzes 

trace  P  <£  P’  (2.4) 

over  all  possible  designs,  where  C4  is  a  generalized  inverse  of  Cj. 

Experience  has  shown  that  this  minimization  is  a  difficult  task.  As 
in  other  cases  of  exact  design  theory,  it  is  highly  unlikely  that  we  can 
obtain  one  method  which  is  capable  of  producing  A-optiaal  designs  for 
arbitrary  values  of  v,  b  and  k.  Recently  several  families  of  A-optiaal 
designs  have  been  discovered. 
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At  this  point  it  1*  useful  to  recall  a  celebrated  result.  If  there 
was  no  control  and  if  ve  were  interested  in  comparing  v  teat  treataents 
aaong  theaselves  than  a  BIB  design  in  the  v  test  treataents  would  be 
A-optlaal.  Unfortunately,  with  the  presence  of  the  control  and  for  the  set 
of  contrasts  of  interest  a  BIB  design  is  alaost  never  an  A-optiaal  design. 
However,  we  can  soaetlaes  utilise  BIB  designs  in  the  teat  treataents  to 
construct  an  A-optiaal  design  for  our  problea.  We  shall  give  soae  faallles 
of  such  designs  below.  For  convenience,  we  introduce  the  notation  ABIB 
(v,b,k-t;t)  to  denote  a  BIB  design  in  the  v  test  treataents  in  b  blocks  of 
else  k-t  each  augaented  by  t  replications  of  the  control  in  each  block. 


Fanlly  1.  An  ABIB(v,b,k-l ;1)  is  A-optimal  whenever  (k-2)*  ♦  1  i  v  i(k-l)*. 

An  exaaple  of  an  A-optiaal  design  when  v  ■  7,  b  ■  7,  k  ■  4  ii  given 
below,  where  the  coluans  are  the  blocks; 


0  0  0  0  0  0  0 

1  2  3  4  5  6  7 

2  3  4  5  6  7  1 

4  5  6  7  1  2  3 

For  each  (v,k)  satisfying  (k  -  2)*  +  1  5  v  6  (k  -  1)*,  there  are  an 
infinite  number  of  A-optiaal  ABIB(v,b,k  -  1 ; 1 )  designs.  These  results  and 
more  details  are  available  in  Hedayat  and  Majumdar  (1985a). 

Stufken  (1986)  has  generalized  the  preceding  idea  to: 


Faaily  2.  An  ABIB(v,b,k  -  t;t)  is  A-optimal  whenever  (k-t  -  1)*  + 
1  S  t2v  S  (k  -  t)». 


An  tuipl*  of  an  A-optiaal  daaign  when  v  ■  8,  b  ■  28,  k  •  8  i»  given 
below: 

0000000000000000000000000000 

0000000000000000000000000000 

3222222111111111111111111111 

4433333433333222222222222222 

5554444554444544443333333333 

6666555666555665556555444444 

7777766777766777667766766555 

8888887888887888878887887876 

Soaetlaes  we  can  use  two  BIB  designs  to  construct  an  A-optiaal  design 
for  our  problea.  Ue  give  below  one  such  f sally,  which  is  taken  froa  Cheng, 
Majuadar,  Stufken  and  Ture  (1986): 

Faally  3.  For  v  •  a*  -  1,  b  ■  y(e  ♦  2)(e*  -  1)  and  k  •  «,  the  union  of  an 
ABIB(v,y(a  +  l)(e2  -  l),a  -  1 ;  1 )  and  a  BIB  design  in  all  the  v  ♦  1  treat- 
aents ,  teat  treataenta  and  control,  in  y«(a  ♦  1)  blocks  of  size  k  each  is 
A-optlaal  whenever  a  is  a  prlae  power,  and  r  is  any  integer. 

An  exaaple  when  a  «  3,  y  •  1,  v  »  8,  b  •  40,  k  ■  3  is: 

00000000000000000000 

11111112222223333344 

23456783456784567856 

00000000147123123123 
44555667258456564645 
78678788360780078807 
Stufken  (1986b)  has  soae  aore  faailies  of  A-optiaal  designs. 


To  establish  tha  optimality  of  theaa  families,  the  starting  point  is  a 
reault  due  to  Majumdar  and  Note  (1983),  which  we  now  proceed  to  state. 
There  are  many  parameter  combinations  (v,b,k)  which  do  not  belong  to  any  of 
the  three  families  but  for  which  Majumdar  and  tiotz's  result  could  still  be 
used  to  get  an  optimal  design.  A  complete  list  of  all  designs  available  by 
this  result,  when  2  £  k  £  8,  k  £  v  £  30,  v  £  b  £  50  is  given  in  Hedayat  and 
Majumdar  (1984).  Before  stating  the  result  we  need  some  definitions. 

Definition  2.1.  d  is  a  Balanced  Treatment  Incomplete  Block  (BTIB)  design  if 

1-do  l  "  «  •  *  *  Ido  v» 

ld»  ^dv-i , v 

b 

where  Xdjj  “  ndipndjp*  This  definition  is  due  to  Bechhofer  and 

Tamhane  (1981). 

Definition  2.2.  For  integers  t  e  (0, 1 , . . . ,k-l)  and  s  t  {0,1 , . . . ,b-l) ,  d  is 
a  BTIB( v,b,k; t ,s)  if  it  is  a  BTIB  design  with  the  additional  property  that 

ndij  *  {0*1) »  1  •  1*  •«•*  v;  J  ■  1,  ...,  b, 

ndo !-.••-  ndo s  *  *  +  1» 

ndo  ,s+i  “  *  ndob  " 

A  BTIB( v,b,k;t,s)  is  called  a  Rectangular  (R-)  type  design  when  a  »  0,  and 
a  Step  (S-)  type  design  when  s  >  0.  The  layout  of  these  designs  can  be 
pictured  as  follows,  with  columns  as  blocks,  in  each  of  the  two  cases 


nents  only 


We  shall  also  naad  the  following  notations: 

A  •  {(x,z):x»0, .. . , (k/2]-l ;*-0, . . . #b  with  z>0  when  x«0)  (2.5) 

a  ■  (v-1)* ,c«bvk(k-l),p-v(k-l)+k 

A(x,z)  ■  {c-p(bx+z)+bx*+2xz+z}/vk 

B(x,z)  -  {k(bx+z)-(bx*+2xz+z)}/vk 

g(x,z)  -  a/A(x,s)  ♦  1/B{x,z). 

Now  we  are  ready  to  state  the  result  of  Majuadar  and  Notz  (1983). 

Theorem  2.1.  Let  v,  b,  k  be  Integers  with  k  1  v.  A  BTIB(v,b,k;t , s)  Is 
A-optimal  In  the  class  of  all  designs  if 

g(t,s)  -  Min  {g(x,z):(xtz)  e  A}. 

Hedayat  and  Majumdar  (1984)  have  devised  an  algorithm  for  obtaining 
A-optimal  designs  based  on  Theorea  2.1.  Jacroux  (1986)  has  generalized 
this  algorithm.  His  algoritha  is  often  capable  of  producing  A-optiaal 
designs  which  are  not  necessarily  BTIB  in  their  structure.  In  particular, 
the  algoritha  given  by  Jacroux  (1986)  often  produces  A-optimal  group 
divisible  treataent  designs  (GOTD's). 

Definition  2.3.  d  is  a  GDTD  with  parameters  n,n,\9,\l  and  \t  if  the 
treatments  1,  ....  v  can  be  divided  into  m  groups  Vt  ,  . ..,  VB  of  size  n 
such  that 

(i)  Idol  "  ^or  *  “  1?  •  ••.  v  and  for  some  constant  X#  , 

(ii)  if  i,j  e  Vp, i#J ,  Xjij  -  Xj  for  some  constant  Xt , 

(iii)  if  i  e  Vp,J  e  Vq,  p#q,  Xjjj  ■  X2  for  some  constant  X2 . 
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Definition  2.4.  For  Integers  t  e  {0,1 , . . . ,k-l}  and  a  e  {0,1 . b-1) ,  d  Is 

•  GDTD(v,b,k;t ,a)  if  it  is  a  ODTD  with  the  additional  property  that 

ndiJ  e  f0,l},  i  ■  l,...,v,J  •  l,...,b, 

|  "  •  •  •  *  ndos  **  t  ♦  1 , 
ndo,s+i  ***••"  ndob  * 

A  GDTD(v,b,k;t,s)  is  called  an  R-type  design  when  s  ■  0  and  an  S-type 
design  when  s>0. 

To  state  Jacroux's  generalization  of  Theorea  2.1,  we  will  let  a,c,p 
A(x,z),  B(x,z),  and  g(x,z)  be  as  defined  in  (2.5)  and  also  introduce  the 
following  notations: 

A  ■  {(x,z):  x  -  0,...,k-2;  z  »  0,...,b  with  z>0  when  x  •  0}  (2.6) 

R(x,z)  •  [(bk  -  bx  -  z)/vj. 

X(x,z)  ■  [ {bk-bx-z)(k-l)  -  vkB(x,z)}/v(v-l) ) 

C(x,z)  -  {bk-bx-z-vR(x,z)){R(x,z)(k-l)}*/k* 

♦{(v-bk+bx+z+vR(x,z)HR(x,z)(k-l)}*/k*+{A(x,z)-v(v-l)l(x,z)HX(x,z) 

♦1}* /k*+{v(v-l)-A(x,z)+v(v-l)X(x,z)}X* (x,z)/k* . 

» (x, z)  ■  {C(x,z)  -  B*(x,z)  -  A* (x,z)/(v-l)}i 
h(x,z)  •  a/(A(x,z)-2/k)  ♦  1/B  (x,z). 

/g(x,z)  if  B(x,z)  >  {A(x,zM(v-l)/(v-2))$P(x,z)}/(v-l), 
a(x,z)  -  <l/B(x,z)+(v-2)(v-l)/{A(x,z)-((v-l)/(v-2))ip(x,z)) 
^♦(v-l)/{A(x,z)+((v-l)(v-2))ip(x,z)},  otherwise 

n(x,z)  ■  ain  {h(x,z),  a(x,z)}. 

Theorea  2.2.  Let  v,  b,  k  be  Integers  with  k  6  v.  A  BTIB(v,b,k;t ,s)  or 
a  GDTD(  v , b , k;  t ,  i )  having  a  •  2,  n  ■  v/2  and  X2  -  X,  ♦  1  is  A-optiaal  in 
the  class  of  all  designs  if 

n(t,s)  ■  ain  {n(x,z)  :  (x,z))  e  A}. 
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We  not*  that  aany  BTIB(v,b,k;t , •)  designs  not  satisfying  the  conditions 

of  Thaoraa  2.1  can  ba  shown  to  satisfy  the  conditions  of  Theorea  2.2.  In 

addition,  Theorea  2.2  can  be  used  to  establish  the  A-optiaality  of 

GDTD(v,b,k;f  ,!)'*  having  a  «  2,  n  •  v/2  and  \t  ■  4-1.  Using  soae 

aore  elaborate  coaputational  techniques,  Jacroux  (1986)  ha*  also  developed 

soae  sufficient  conditions  for  GDTD(v,b,k;t  ,s)'s  having  lt  *  1,  4  1  or 

a  ■  v/2,  n  ■  2  and  lt  ■  1,  ■  1  to  be  A-optlaal  aaong  all  designs  with 

paraaeters  v,  b,  and  k.  One  exaaple  of  an  A-optlaal  GOTO  is  that 

GDTD(9,9,4;1,0)  d*  having  a  •  3,  n  •  3,  \j  ■  0  and  lt  •  1  given  below: 

000000000 
1  1  1  2  2  2  3  3  3 

456456456 
789897978 


2.2.  A-optiaal  designs  for  2-way  eliainatlon  of  heterogeneity.  Our 
statistical  setup  consists  of  bk  experimental  units  srranged  in  a  k  x  b 
array,  and  the  model  of  response  under  design  d  is: 

ydijjt  *  P  ♦  *i  ♦  ®J  ♦  PJL  ♦  eijjl  <2-7> 

(i  ■  0,  1,  ....  v;  J  »  1,  ...,  b;  %  -  1,  ...,  k)  if  treatment  1  is  spplled 
to  the  experimental  unit  in  cell  (l,j). 

Let 

ndij  “  number  of  times  treatment  i  occurs  in  column  j , 

mdi£  *  number  of  times  treatment  i  occurs  in  row  1, 

b 

rdi  *  ndi j  » 

Nd  *  <ndiJ7*  •  v  x  b  matrix, 

*d  *  (“dip*  a  v  x  k  matrix, 

P  is  the  v  x  (v  ♦  1)  matrix  defined  in  subsection  2.1, 
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Cj(j)«  (  rdt  •  rdi  '  ’  *  *  •  rdv^ ”k_1  ♦  (bk)"‘rjrj  . 

Than  an  A-opfciaal  design  ilaiiltai  traea  PC^^P'.  Ha  shall  now  high¬ 
light  iom  of  tha  rasulta  f ro«  recant  literature. 

Family  1.  Let  p  be  an  integer  and  v  ■  p* ,  b  ■  k  ■  p*  ♦  p.  Abxb  array 
In  which  each  taat  treatment  appears  once  in  each  row  and  in  each  column 
and  the  control  appear*  p  time*  in  each  row  and  in  each  column  la  A- 
optlmal . 

One  eaay  way  to  construct  members  of  this  family  is  to  start  with  a 
Latin  square  of  order  p*  ♦  p  and  change  symbol*  p*  ♦  1,  . . . ,  p*  ♦  p  to  0 
(control).  He  illustrate  this  in  the  following  example  with  v  ■  4,  b  -  k  ■ 
6. 
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This  and  some  more  general  results  are  available  in  Hots  (1985). 

This  has  been  generalized  by  Majumdar  (1986). 

Family  2.  Let  p,  a  and  7  be  integers,  v  ■  p*  ,  k  •  o(p*  ♦  p)  and  b  • 

7<p2  +  p).  A  k  x  b  array  in  which  each  test  treatment  appears  o  times  in 
each  column  and  y  times  in  each  row,  and  the  control  appears  op  times  in 
each  column  and  yp  times  in  each  row  is  A-optimal. 


-15- 


One  way  to  construct  atabari  of  this  faaily  is  to  fora  tho  array 
(L^j),  i  «  1,  .  .  .,  a;  j  »  1,  .  .  .  y 
whoro  aach  L^j  la  a  atabar  of  Faailly  1. 

Faaily  3.  A  k  x  b  array  la  A-optlaal  if 

(1)  it  la  an  A-optlaal  block  doaign  for  1-vay  ellalnation  of 
hoteroganoity  with  columns  aa  blocka,  and 
(11)  tha  total  nuabar  of  raplicationa  for  aach  treatment,  taat 
treataant  or  control,  la  divided  equally  aaong  tha  k  rowa. 

Thla  haa  baen  given  by  Jacroux  (1984c).  The  following  ia  an  exaapla 
whan  v  ■  9,  b  ■  24,  k  ■  3: 

041025783096008079162345 

105802005204466907213879 

310140220730950680976458 

3.  MV-optiaal  Designs.  In  thla  aaction  we  give  a  nuabar  of  raaulta 
concerning  tha  MV-optlaality  of  dealgna  for  comparing  v  taat  treataenta 
with  a  control  for  tha  0-way,  1-way  and  2-way  ellalnation  of  heterogeneity. 
The  notation  Introduced  in  Section  2  ia  alao  uaed  throughout  thla  section. 

3.0.  MV-optlaal  dealgna  for  0-way  aliaination  of  heterogeneity.  Our 
atatlatlcal  setup  is  tha  sane  as  given  in  Subsection  2.0.  An  MV-optiaal 
design  minimizes 

aax  (l/rj,  +  l/r^)  (3.1) 

iSiiv 

subject  to  the  restriction  r<j#  ♦  r{j1  ♦  . . .  ■♦■  r<jv  ■  n.  Using  elementary 
mathematics,  it  is  easily  seen  that  an  MV-optimal  design  d*  has 
rd*i  "  *«+lfor  i  -  1,  ....  v  and  r^*#  ■  n  -  + 


where  f|  la  given  in  (2.3).  We  note  that  for  e  fixed  value  of  n,  the 
A-and  the  MV-optimality  criteria  may  select  substantially  different  optimal 
designs  from  those  available.  For  example,  when  n  ■  30  and  v  ■  15,  an 
A-optlmal  design  d  will  have  rjj*  *  5  and  rj^  •  1  or  2  for  i  ■  1,  .... 
15  whereas  the  MV-optlswl  design  d*  will  have  r^  ■  15  and  ■  1  for  1 
*  1 ,  . . . ,  15. 

3.1.  MV-optimal  designs  for  the  1-way  elimination  of  heterogeneity.  Our 
statistical  set-up  is  the  same  as  given  in  Subsection  2.1.  We  note  that 
any  design  which  is  A-optlmal  among  all  designs  having  parameters  v,  b  and 
k  and  which  estimates  all  contrasts  of  the  form  t±  -  t,  with  the  same 
variance  will  also  be  MV-optimal.  Thus  we  see  that  all  designs  given  in 
Subsection  2.1  as  being  A-optlmal  are  also  MV-optimal  since  all  GDTD 
( v,b, k; T, 5) ' s  estimate  contrasts  of  the  form  tj  -  tt  with  the  same 
variance.  However,  Jacroux  (1986)  has  developed  some  additional  sufficient 
conditions  which  can  be  used  to  establish  the  MV-optlsiallty  of  various 
GDTD( v,b,k;x, z) ' s  which  cannot  be  proven  to  be  A-optimal  using  any  known 
results.  As  an  example  of  the  types  of  results  which  can  be  proven  for 
MV-optlmality,  we  have  the  following. 

Theorem  3.1.  Let  d*  be  a  BTIB(v,b,k;t,s)  where 

n(t,s)  •  min  {n(x,z) : (x,z)  c  A  ,  (bk  -  bx  -  z)/v  is  an  integer) 
and  for  positive  integers  p  and  q,  let  B  (p,q)  denote  the  smallest  value 
of  y  such  that 

vhtrt  t  *  it  the  lust  squirts  estlaiti  of  t.  obtained  under  d*. 
a  i  i 
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Alio  lat  R(x,s)  •  f(bk  -  bx  -  t)/v).  If  for  any  (x,x)  c  A  inch  that 
n(x,x)  <  n(t,a),  It  hold*  that 

vkB(x,s)  <  <bk-bx-*-vR(xfa))B(vkB(x,i) ,(R(x,z)+l)(k-l)> 
+(v-bk+bx+*+vR(xf *))B(vkB(x,r) ,R(x,*)(k-l)) , 
than  d*  ia  MV-optimal  aaong  all  daalgna. 

Ualng  aoaa  aora  coaplax  coaputatlonal  tachniquaa,  Jacroux  (1986)  haa 
obtalnad  aoaa  further  raaulta  aiailar  to  Thaoraa  3.1  which  can  ba  uaad  to 
eatabliah  tha  MV-optlaallty  of  varioua  GDTD(v,b,k;t,a) 'a  having  0  6 
*2  “  11  or  a  ■  v/2  n  -  2  and  ■  1,  -  1  vhoaa  A-optlaallty  raaalna 

unknown . 

For  axaapla  whan  v  •  6,  b  ■  11  and  k  ■  3,  aa  wall  aa  whan  v  ■  6,  b  ■ 
16  and  k  •  4,  an  A-optlaal  daalgn  ia  unknown.  However,  wa  ara  able  to  glva 
daalgna  whoaa  MV-optlaallty  can  ba  aatabllahed  ualng  raaulta  auch  aa 
Thaoraa  3.1.  Thaaa  ara  axhlbltad  below: 

Exaapla  3.1.  v  ■  6,  b  •  11  and  k  ■  3. 

0000000001  4 

11122233325 
45645645646 


It  la  interesting  to  note  that  the  design  in  Example  3.1  is  an 
S-type  BTIB  design,  while  the  design  given  in  Example  3.2  is  a  GDTD. 

3.2.  MV-optlmal  designs  for  the  2-way  elininatlon  of  heterogeneity. 
Again  our  statistical  set-up  is  the  same  as  that  given  in  Subsection  2.2. 
Using  arguments  similar  to  those  used  in  Subsection  3.1,  we  see  that  all  of 
the  A-optimal  row-column  designs  which  estimate  treatment  contrasts  tj-tf 
with  the  sane  variance  will  also  be  MV-optlmal.  Thus  all  the  A-optinal 
row-column  designs  listed  in  Subsection  2.2  are  also  MV-optlmal.  In 
addition,  a  k  x  b  array  is  MV-optinal  if 

(1)  it  is  an  MV-optlmal  block  design  for  1-way  elimination 
of  heterogeneity  with  columns  as  blocks,  and 
(11)  the  total  number  of  replications  for  each  treatment,  test 
treatment  or  control,  ia  divided  equally  among  the  k  rows. 

Example  3.3.  For  v  ■  6,  b  ■  16  and  k  ■  4,  the  row-column  design  given  by 

0645031  205620314 

1  0  64205310435062 

2103460632051405 
3520154046102630 

is  MV-optlmal. 

4.  Model  Robust  Optimal  Designs.  There  are  circumstances  in  which  the 
experimenter  is  not  sure  whether  to  fit  a  1-way  or  a  2-way  elimination  of 
heterogeneity  model  to  the  data.  In  such  situations  it  would  be  highly 
desirable  to  obtain  a  design  which  is  A-  or  MV-optlmal  under  each  of  these 


models.  Rtdtyat  and  Majumdar  (1986b)  atudiad  thla  aapact  of  tha  problaa 
and  gava  aoaa  faniliaa  of  aodal  robuat  dealgna.  Tha  faailiaa  were  con- 
atruetad  using  tha  Euclidean  plana,  the  Projective  plane  and  soae  other 
geoaatrical  structures.  The  exact  description  of  the  faallies  are  soaevhat 
involved;  soae  typical  examples  are  given  below. 


Example  4.1.  Let  v  •  4,  k  ■  3  and  b  ■  6.  The  following  design  is  A-  and 
MV-optimal  for  both  1-  and  2-way  elimination  of  heterogeneity  models: 


In  fact,  this  design  is  A-  and  MV-optimal  for  the  O-way  elimination  of 
heterogeneity  model  as  well. 


Example  4.2.  Let  v  ■  7,  k  ■  4  and  b  •  28.  The  following  design  is  A-and 
MV-optimal  for  both  1-  and  2-way  elimination  of  heterogeneity  smdels: 


0000000  1234567  1234567  1234567 
1234567  0000000  2345671  2345671 
2345671  2345671  0000000  4567123 
4567123  4567123  4567123  0000000 


Before  closing  this  section  we  would  like  to  mention  that  the  designs 
in  families  1  and  2  in  subsection  2.2  are  A-  and  MV-optimal  under  0-way, 
1-way  and  2-way  elimination  of  heterogeneity  models,  while  the  designs  in 
family  3  are  A-  and  MV-optimal  at  least  under  1-way  and  2-way  elimination 
of  heterogeneity  models. 
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5.  OthT  Efficient  D>»ign».  Evan  though,  for  each  set  of  v  teat  treatment 
there  ie  an  A-  or  MV-optimal  design  for  a  0-,  1-  or  2-way  elimination  of 
heterogeneity  model,  the  task  of  finding  this  design  can  be  very  difficult 
indeed.  For  situations  where  an  A-  or  MV-optinal  design  is  unknown,  there 
are  several  alternative  ways  of  planning  an  experiment.  Here  are  aoae 
possibilitlea. 

5.1.  Liait  the  class  of  competing  designs  to  a  "reasonably  rich"  aubclaaa, 
so  that  an  A-  or  MV-optimal  design  within  this  subclass  can  be  constructed. 
For  example,  under  a  1-way  ellainatlon  of  heterogeneity  model  the  BTIB 
designs  fora  such  a  subclass.  For  v  ■  3,  b  ■  15,  k  »  2  the  A-  or  MV- 
optimal  design  is  not  available  in  the  literature,  but  a  design  which  is  A- 

I 

and  MV-optimal  among  all  BTIB  designs  is  given  by 

'  000000000112 
111222333233 

This  approach  has  been  studied  in  Hedayat  and  Majumdar  (1984)  under  the 

A-optlaallty  criterion,  and  some  series  of  such  designs  have  been 

cataloged.  Numerical  evidence  indicates  that  optimal  designs  obtained  in 

this  fashion  are  highly  efficient  in  the  class  of  all  designs.  There  are, 

however,  Isolated  Instances  where  they  perform  poorly.  A  similar  study  for 

the  MV-optlmallty  criterion  has  been  carried  out  by  Jacroux  (1985). 

5.2.  Search  for  an  approximately  A-  or  MV-optimal  design.  This  can  be 
carried  out  in  two  steps.  First  compute 

(i)  g(t,s)  -  Min{g(x,z) : (x,z)  c  A) 

or  (5.1) 

(ii)  g(t,s)/v  »  Min{g(x,z) :(x,z)  e  A}/v 
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where  the  function  g(x,s)  has  boon  defined  in  Subsection  2.1.  Those  give 

lower  bounds  to  the  values  of  the  A-  and  MV-criterion,  respectively.  The 

second  step  consists  of  guessing  a  good  design  d  based  on  available 

theory.  Coapute  the  corresponding  value  of  the  expression  (2.4)  or 

aax  var  (t-H-tH.i  for  this  design  and  coapare  with  the  appropriate  alnlaua 
ISiSv 

value  given  in  (5.1).  If  the  coaparison  la  poor  in  the  opinion  of  the 

experiaenter,  then  he  should  aodlfy  his  guess  and  try  again. 

Let  us  deaonstrate  this  approach  by  an  exaaple.  Let  v  •  21,  b  »  30 

and  fc  -  9.  Here  the  ainiaua  given  by  (5.1(1))  is  2.589  and  that  given  by 

(5.1(11))  is  .1233.  Our  experience  shows  that  BIB  designs  in  the  test 

treataents  augaented  by  one  or  aore  replications  of  the  control  in  each 

block  are  often  highly  efficient,  as  seen  in  faailies  1  and  2  of  Subsection 

2.1.  In  our  case  we  can  try  a  design,  d ,  which  is  an  ABIB(21 , 30, 7;2) . 

For  this  design  the  value  of  (2.4)  is  2.618  and  aax  var(tjj  -  tj.^  “  *1247, 

ISiSv 

giving  an  efficiency  of  at  least  98.87%  for  both  the  A-  and  MV-optiaality 
criteria.  So  this  is  Indeed  a  highly  efficient  design.  This  approach  of 
approxiaating  an  A-and  MV-optinal  design  by  an  augaented  BIB  design  has 
been  studied  by  Stufken  (1986c). 

Another  aethod  of  tracking  down  a  good  approxiaatlon  has  been  given 
in  Cheng,  Majuadar,  Stufken  and  Ture  (1986).  It  consists  of  first 
deteraining  the  point  (t,s)  which  minimizes  the  function  g(x,z)  given  in 
(5.1).  In  case  a  BTIB(  v ,  b  ,  k ;  t ,  s )  exists,  it  is  both  A-  and  MV-optiaal. 

If  it  does  not,  then  at  least  one  of  the  following  two  designs  is  expected 
to  be  a  good  approximation: 

(i)  A  design  with  the  same  number  (bt+s)  of  replications  of  the  con¬ 
trol  as  a  BTIB(v,b,k; t , s)  and  which  is  "combinatorially  close" 
to  a  BTIB  design. 


(11)  A  BTZB  design  with  the  number  of  replications  of  the  control 
"close"  to  bt  ♦  e. 

We  demonstrate  the  idea  by  an  exaaple  when  v  ■  5,  b  ■  7  and  k  ■  4. 
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Consider  the  following  two  designs: 

0 

0  0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1  1 

1 

1  2 

2 

0 

0 

1 

1 

2 

2 

2 

d»  S  2 

2  3 

3 

4  3 

3  * 

1 

3 

3 

4 

3 

4 

3 

4 

5  4 

5 

5  4 

5 

2 

4 

5 

5 

5 

5 

4 

Here  bt  ♦  s  ■  7,  d,  is  a  non-BTIB  design  with  7  replications  of  the 
control,  while  dt  is  a  BTIB  design  with  8  replications  of  the  control. 
The  value  for  expression  (2.4)  for  d,  is  2.058  giving  it  an  efficiency 
of  99. 2%.  The  value  for  expression  (2.4)  for  dt  is  2.143  giving  it  an 
efficiency  of  95. 2Z.  Thus  both  of  these  designs  are  highly  efficient 
under  the  A-crlterion,  dt  being  the  better  of  the  two. 

Finally,  with  the  availability  of  today's  high-speed  computers  and 
supercomputers,  one  can  find  an  A-  or  MV-optimal  design  by  a  complete 
search  among  all  possible  designs  if  the  parameters  v,  b  and  k  are  not  too 
large. 


6.  A-  and  MV-optimal  Designs  for  Two  or  More  Controls.  So  far  we  have 
been  discussing  optimal  designs  for  comparing  test  treatments  with  one 
control.  There  are  circumstances  when  the  test  treatments  have  to  be 
compared  with  two  or  more  controls.  Suppose  we  denote  by  S  the  set  of  all 
controls  and  by  T  the  set  of  all  test  treatments.  Then,  an  A-optimal 
design  is  the  one  which  minimizes 


> 
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gfS  hSl  var<?d*  *  *dh>. 

and  an  MV-optiaal  design  is  ona  which  minimizes 

_  var(tdg  -  tdh) 
gcS.heT  * 

aaong  all  designs  under  a  O-way,  1-way  or  2-way  elimination  of  hetero¬ 
geneity  aodel.  Najuadar  (1986)  has  studied  this  problea,  and  has  identi¬ 
fied  designs  which  are  A-  and  MV-optiaal  in  various  settings.  Here  we 
give  an  example  of  an  A-  and  MV-optiaal  design  for  one  way  elimination  of 
heterogeneity  and  an  example  of  a  design  which  is  optimal  for  each  of  0- , 
1-  and  2-way  elimination  of  heterogeneity. 

Example  6.1.  Suppose  4  test  treatments  are  to  be  compared  with  3  controls 

in  30  blocks  of  size  3  each  under  the  1-way  elimination  of  heterogeneity 

model.  Denoting  the  test  treataients  by  A,B,C,D  and  the  controls  by  1,  2, 

and  3  the  A-  and  MV-optimal  design  is: 

112112112112111111 

233233233233AAABBC 

AAABBBCCCDDDBCDCDD 

222222333333 

AAABBCAAABBC 

BCDCDDBCDCDD 

Example  6.2.  Suppose  8  test  treatments  are  to  be  compared  with  2  controls 
in  a  12  x  12  array  under  a  2-way  elimination  of  heterogeneity  model. 
Denoting  the  test  treatments  by  A,B,C,D,E,F,G,H,  and  the  controls  by  1,2, 
the  A-  and  MV-optimal  design  is  given  by 
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A  B 
B  C 

C  D 
D  E 
E  F 
F  G 
G  H 
H  1 
1  1 
1  2 
2  2 
2  A 


C  0 
D  E 
E  F 
F  G 
G  H 
H  1 
1  1 
1  2 
2  2 
2  A 
A  B 
B  C 


E  F 
F  G 

G  H 
H  1 
1  1 
1  2 
2  2 
2  A 
A  B 
B  C 
C  D 
D  E 


G  H 
H  1 
1  1 
1  2 
2  2 
2  A 
A  B 
B  C 
C  0 
D  E 
E  F 
F  G 


1  1 
1  2 
2  2 
2  A 
A  B 
B  C 
C  D 
D  E 
E  F 
F  G 
G  H 
H  1 


2  2 
2  A 
A  B 
B  C 
C  D 
D  E 
E  F 
F  G 
G  H 
H  1 
1  1 
1  2 


This  design  is  model  robust  in  the  sense  of  being  A-  and  MV-optlmal  for  0— 
and  1-vay  elimination  of  heterogeneity  models  as  well. 


7.  Discussion.  Cox  (1958,  p.238)  advocated  augmenting  a  BIB  design  in 
test  treatments  with  one  or  more  replications  of  controls  in  each  block  as 
a  means  of  getting  good  designs.  He  neither  formally  mathematized  the 
problem  nor  gave  any  justification  for  his  suggestion.  However,  based  on 
what  has  been  developed  during  the  past  several  years,  we  know  that  this  is 
an  excellent  method  of  getting  efficient  designs  in  many  cases.  Pearce 
(1960)  gave  a  solution  for  A-optimal  designs  for  the  0-way  elimination  of 
heterogeneity  model  which  is  applicable  when  v  is  a  square,  and  proposed  a 
class  of  designs  for  comparing  test  treatments  with  a  control  and  gave 
their  analysis  for  the  1-way  elimination  of  heterogeneity  model.  Freeman 
(1975)  studied  some  designs  for  comparing  two  sets  of  treatments  for  the 
2-way  elimination  of  heterogeneity  model.  Pesek  (1974)  compared  a  BIB 
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design  with  sn  augmented  BIB  design,  as  suggested  by  Cox  (1958),  in 
estimating  control-test  treatment  contrasts  and  noticed  that  the  latter  was 
more  efficient.  Das  (1958)  has  also  looked  at  augmented  BIB  designs. 

Bechhofer  and  Tamhane  (1981)  were  the  first  to  study  the  problem  of 
obtaining  optimal  block  designs.  However  their  optimality  consideration 
was  neither  A-  nor  MV-optimality,  but  for  the  problem  of  obtaining  optimal 
simultaneous  confidence  intervals  under  a  1-way  elimination  of  hetero¬ 
geneity  model.  Their  discoveries  led  to  the  concept  of  BTIB  designs;  Notz 
and  Tamhane  (1983)  studied  their  construction. 

Constantine  (1983)  showed  that  a  BIB  design  in  test  treatments 
augmented  by  a  replication  of  the  control  in  each  block  is  A-optimal  in  the 
class  of  designs  with  exactly  one  replication  of  the  control  in  each  block. 
Jacroux  (1984a)  showed  that  Constantine's  conclusion  remains  valid  even 
when  the  BIB  designs  are  replaced  by  some  group  divisible  designs. 

Majumdar  and  Notz  (1983)  gave  a  method  of  obtaining  A-  and  MV-optimal 
designs  among  all  designs  for  the  1-way  elimination  of  heterogeneity  model. 
Hedayat  and  Majumdar  (1984)  gave  an  algorithm  and  a  catalog  of  A-  and 
MV-optimal  designs  and  studied  approximations.  Ture  (1982)  also  studied 
A-optlmal  designs  and  their  approximations  and  construction.  Hedayat  and 
Majumdar  (1985a)  gave  families  of  A-  and  MV-optimal  designs.  Notz  (1985) 
studied  optimal  designs  for  the  2-vay  elimination  of  heterogeneity  model. 
Majumdar  (1986)  considered  the  problem  of  finding  optimal  designs  for 
comparing  the  test  treatments  with  two  or  more  controls. 

Jacroux  (1984b,  1985)  gave  new  methods  for  obtaining  MV-optimal 
designs  under  1-way  elimination  of  heterogeneity  models,  gave  catalogs  and 
studied  approximations.  Jacroux  (1984c)  studied  optimal  designs  for  2-way 
elimination  of  heterogeneity  models,  utilizing  techniques  of  Hall  (1935) 
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and  Agarwal  (1966).  Hedayat  and  Majuadar  (1985b)  atudled  designs  simul¬ 
taneously  optlaal  under  both  1-  and  2-way  ellainatlon  of  heterogeneity 
aodels.  Jacroux  (1986)  generalized  the  Hedayat  and  Majuadar  (1984) 
algoritha  for  finding  A-optlaal  designs.  Cheng,  Majuadar,  Stufken  and  Ture 
(1986)  gave  new  faailies  of  A-  and  MV-optiaal  designs  and  soae  approxima¬ 
tions  for  1-way  ellainatlon  of  heterogeneity  aodels.  Stufken  (1986a, b,c) 
studied  A-  and  MV-optlaal  designs  for  1-way  ellainatlon  of  heterogeneity 
aodels,  gave  faailies  and  studied  approxiaations. 

There  are  aany  other  design  settings  in  which  it  would  be  useful  to 
identify  optlaal  designs  for  comparing  test  treatments  with  controls.  One 
such  setting  is  that  of  repeated  measurements  designs.  Soae  aspects  of 
optiaallty  and  construction  of  designs  in  this  area  has  been  investigated 
by  Pigeon  (1984). 

Glovagnoll  and  Wynn  (1985)  studied  A-optlmallty  of  designs  for  1-way 
elimination  of  heterogeneity  aodels  set  in  the  context  of  approxiaate 
theory,  i.e.  with  an  infinite  nuaber  of  observations.  Spurrier  and  Edwards 
(1986)  did  a  siailar  study  for  optimal  designs  for  finding  sisiultaneous 
confidence  intervals. 

It  seems  appropriate  to  make  a  comment  on  randomization.  In  running 
optimal  designs  we  often  have  to  follow  a  well  structured  pattern.  This 
does  not,  however,  mean  that  there  will  be  no  room  for  randomization.  The 
labelling  of  the  treatments,  experimental  units  under  a  O-way  elimination 
of  heterogeneity  model,  blocks  under  a  1-way  elimination  of  heterogeneity 
model  and  rows  and  columns  under  a  2-way  elimination  of  heterogeneity  model 
can  be  randomized. 
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